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The present spatial distribution of galaxies in the Universe is non-Gaussian, with 40% skew- 
ness in 50 h~ Mpc spheres, and remarkably little is known about the information encoded in it 
about cosmological parameters beyond the power spectrum. In this work we present an attempt 
to bridge this gap by studying the bispectrum, paying particular attention to a joint analysis with 
the power spectrum and their combination with CMB data. We address the covariance properties 
of the power spectrum and bispectrum including the effects of beat coupling that lead to interest- 
ing cross-correlations, and discuss how baryon acoustic oscillations break degeneracies. We show 
that the bispectrum has significant information on cosmological parameters well beyond its power 
in constraining galaxy bias, and when combined with the power spectrum is more complementary 
than combining power spectra of different samples of galaxies, since non-Gaussianity provides a 
somewhat different direction in parameter space. In the framework of flat cosmological models we 
show that most of the improvement of adding bispectrum information corresponds to parameters 
related to the amplitude and effective spectral index of perturbations, which can be improved by 
almost a factor of two. Moreover, we demonstrate that the expected statistical uncertainties in erg 
of a few percent are robust to relaxing the dark energy beyond a cosmological constant. 
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I. INTRODUCTION 

Several recent studies have stressed the importance of 
combining different observations to constrain cosmologi- 
cal parameters. A clear example is provided by the anal- 
ysis of the galaxy power spectrum in the Sloan Digital 
Sky Survey (SDSS) jlj, and in the 2dF Galaxy Redshift 
Survey 0|, which have shown the central role played by 
the information contained in the large-scale galaxy dis- 
tribution to break the degeneracies still present in the 
cosmic microwave background (CMB) data despite the 
precision of the WMAP satellite observations |3, U ■ 

One of the main challenges in extracting cosmological 
information from galaxy clustering is knowing how good 
tracers of the underlying mass distribution galaxies are. 
This is often bypassed altogether, for example in P, Q 
only infomation on the shape of the galaxy power spec- 
trum was used, since its amplitude is degenerate with 
the linear bias parameter relating galaxy to dark matter 
fluctuations at large scales. 

The determination of galaxy bias has been, 
among the main reasons of interest in the galaxy 
order statistics in general 0, S S IS HOi E 
04, ^51 llq and the bispectrum in particular []J 
123, Izj- At large scales, the dependence on triangle con- 
figuration of the bispectrum generated by gravitational 
instability allows to disentagle the gravitational contribu- 
tion from the bispectrum generated by non-linear biasing 
and ultimately remove the degeneracy between the linear 
bias and the amplitude of the dark matter perturbations. 
In weak gravitational lensing at smaller scales, the bis- 
pectrum can similarly be used to break degeneracies be- 
tween matter content and the amplitude of fluctuations 
and probe dark energy [H El, H El El . 

Observational applications of this method to the 
galaxy distribution in t he p ast involved fixing the cos- 
mological model [la. HM ligl 1201 , thus the information of 




the bispectrum on cosmological parameters has not been 
properly taken advantage of. In this work we study the 
constraining power of the bispectrum as a tool in the de- 
termination of cosmological parameters and the nature of 
primordial fluctuations, going beyond the determination 
of galaxy bias alone. As shown in |60j |. higher-order cor- 
relation functions such as the bispectrum or the trispec- 
trum in galaxy surveys show, when all measurable con- 
figurations are taken into account, a signal-to-noise ratio 
comparable or even exceeding the signal-to-noise of the 
power spectrum at mildly non-linear scales. 

Postponing a detailed discussion of our method to the 
following sections, we show in Fig.^how the power spec- 
trum and the bispectrum measured from the same data 
set compare in constraining cosmological parameters. We 
consider flat cosmological models depending on nine pa- 
rameters: the physical dark matter density ujd = ^dh^, 
the physical baryon density Wf, = Q,i,h?, the dark energy 
density f^A, the amplitude of scalar fluctuations Ag, the 
scalar spectral index n^, the dark energy equation of state 
parameter w, the optical depth to Thomson scattering 
T, plus the linear (6i) and quadratic (62) galaxy bias pa- 
rameters. We also show "derived" parameters such as 
h, the Hubble constant in units of 100km s^^ Mpc~^, 
the baryon density fib and the amplitude of dark matter 
fluctuations at 8/i~^Mpc, a%. 

These constraints are from a hypothetical analysis 
that combines the CMB data from WMAP (first year) 
with measurements in the North part of SDSS by the 
end of the survey in two cases: using the SDSS galaxy 
power spectrum (blue, dashed line) and replacing the 
SDSS power spectrum by the bispectrum (red, contin- 
uous line). Both cases include the covariance between 
different power spectrum bins or bispectrum configura- 
tions (see below for a full discussion). Figure ^ shows 
that when all triangle configurations are included down 



to wavenumber k„ 



0.3 /i Mpc the bispectrum even 
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FIG. 1: An example showing the constraining power of the 
bispectrum compared to the power spectrum. The panels 
show marginalized likelihood functions corresponding to a 
hypothetical joint analysis of WMAP (first year) and SDSS 
North (by the end of the survey) where only the galaxy power 
spectrum is used (blue, dashed line), or only the galaxy bis- 
pectrum is used (red, solid line). Assumes a flat cosmology 
and scales up to k„,.^„ = 0.3/iMpc~^. 



improves the power spectrum results. 

In practice one would like to combine the information 
in the power spectrum and bispectrum, which requires 
a calculation of their covariance properties. This is the 
main subject of the present work. The cross-covariance 
between the power spectrum and bispectrum turns out 
to have some non-trivial properties that help constrain- 
ing cosmological parameters, in a way that is unexpected 
from a separate consideration of the covariance of each 
statistic by itself. 

Although the main properties of the covariance matri- 
ces can be understood analytically, the details of the sur- 
vey under consideration are important in practice, thus 
we compute covariance matrices from mock catalogs de- 
signed to reproduce the geometry of the SDSS survey 
by its completion. In particular, we consider the power 
spectrum and bispectrum of the north hemisphere main 
sample (MS) of galaxies. We also discuss how our con- 
strains change as we include the power spectrum of the 
Luminous Red Galaxies (LRG) sample in the same ge- 
ometry. As an example of what should be expected in 
combining large-scale structure (LSS) with the CMB we 
use the WMAP first year data; after the present work ap- 
peared as a preprint the WMAP 3-year (WMAP3) data 
became available Q . The new data provides a somewhat 
different angle on the constraining power of the bispec- 



trum, not just incremental improvements, for this reason 
we consider separately in the Appendix what happens 
when WMAP3 is added. 

This paper is organized as follows. In section ^ we 
review some basic results regarding the large-scale bis- 
pectrum of the galaxy distribution and discuss the main 
features of the covariance measured from our mock cat- 
alogues in section IIIII In section IIVI we present the 
likelihood functions both for the LSS and CMB corre- 
lations. In scction[3we present results on expected con- 
straints on cosmological parameters, derived assuming 
the WMAP 1-year data (WMAPl). We conclude in sec- 
tion ^T] The Appendix discusses what happens when we 
replace WMAPl by WMAP3. 



II. PREDICTIONS AND MOCK CATALOGS 

We will assume primordial fluctuations to be Gaussian, 
so that every connected higher-order correlation function 
in the dark matter overdensity field S results from gravi- 
tational instability. The dark matter bispectrum B, i.e. 
the Fourier counterpart of the 3-point correlation func- 
tion, is defined as the ensemble average 



('5ki'5k2'^k3> = '5£>(ki23) B{ki,k2,k3), 



(1) 



with (5k the density contrast in Fourier space and ki23 = 
ki + k2 + ks. If the bispectrum can be reliably predicted 
by tree-level perturbation theory (PT), it follows that [23 



B(fci, k2,k3) = 2i^2(ki, k2)PiF2 + eye. 



(2) 



where Pi — P{ki) is the linear power spectrum while the 
kernel F2 reads 

with a; = ki • k2. 

A second source of non-Gaussianity in the galaxy den- 
sity field is given by non-linear galaxy bias. At scales 
much larger than the typical size of virialized structures 
the relation between the galaxy distribution and the un- 
derlyin g da rk matter distribution is expected to be lo- 
cal j29|. I30I l31|. that is, in terms of the respective over- 
densities, 5g{x) = f[S{x)]. For small fluctuations we can 
Taylor-expand and describe such function in terms of few, 
constant, bias parameters |31| 



Sg{x) = bi5 + 



1, 



k 



(4) 



The large-scale galaxy power spectrum will therefore be 
given Pg{k) ~ b\P{k) while the galaxy bispectrum Bg 
will be related to dark matter bispectrum B as 

Bgiki,k2,k3) ~ blB{ki,k2,k3)+blb2iPlP2+CyC.) (5) 

The different behaviour of the first and second terms on 
the right hand side of Eq. jSJ) as a function of triangle 



configuration given by the wavenumbers fci, k2 and k^ 
allows a simultaneous measurement of the linear bias pa- 
rameter 61 and the quadratic bias parameter 62 0, UM ■ 
This becomes obvious when Eq. |(SJ) is rewritten in terms 
of the reduced bispectrum, defined, for the dark matter 
field as Q = B{ki,k2, kz)/{PiP2 + eye.) and analogously 
for the galaxy distribution, so that 



Qgiki,k2,k3) = —Q{ki,k2,k3) + -^ 



bi 



bl 



(6) 



While the first term on the left hand side depends on the 
specific triangle via the F2(ki, k2) kernel, the second just 
amounts to an overall additive constant. 

As mentioned above, the scale-dependence of the bias 
parameter is expected to be very weak at the scales con- 
sidered in the present analysis. This can be shown in the 
framework of the halo model and it can probed, observa- 
tionally, by looking at the dependence of measured values 
of 61 and &2 on the smallest scale, or largest wavenumber 
fc^j^x, included in the analyis. If there is a scatter about 
the deterministic relationship given by Eq. Q), the bis- 
pectrum method recovers the mean relationship between 
galaxy and matter overdensities. This has been shown 
for models with significant scatter (see Fig. 1 in 53] ) and 
galaxies populated using Halo Occupation Distributions 
(HOD) where the scatter is typically not very significant 
at the scales we consider here (see Fig. 6 in oZ\). 

In this work we consider scales up to fc < 0.3 /i Mpc^^, 
for which the validity of Eq. Q is only accurate to about 
20% |5l|, [53 • A more accurate description of the bis- 
pectrum at these scales, particularly in rcdshift space, 
is given by second-order Lagrangian PT (2LPT), |53l| . 
Therefore, we will only use tree-level PT to model devia- 
tions from a fiducial model calculated by using mock cat- 
alogs generated by PTHalos [23 and 2LPT simulations, 
which are similar at these scales, since halos in PTHalos 
are placed in the large-scale 2LPT density field. The 
advantage of using PTHalos is that a biased population 
of galaxies can be chosen by using appropriate prescrip- 
tions for their occupation inside halos. This method is 
therefore necessary for LRG galaxies, which are strongly 
biased tracers, whereas main sample galaxies are close 
enough to being unbiased that the difference between us- 
ing 2LPT and PTHalos is not significant. 

We therefore use the mock catalogs for the main 
sample of galaxies in SDSS generated by using 2LPT 
in |59| , where the following cosmological parameters were 
used: dark matter density Vtd — 0.225, baryon den- 
sity f^b = 0.045, cosmological constant with density 
JIa = 0.73, Hubble constant h = 0.71 and fluctuation 
amplitude a% — 0.82 at the mean redshift of the sur- 
vey of Zmcan — 0.1. As discusscd abovc, we assume 
these galaxies to be unbiased, and have included the de- 
tailed geometry of the expected final angular and radial 
selection functions. The redshift-space density field is 
weighted using the Feldman-Kaiser-Peacock (FKP) pro- 
cedure m il m with Po = 5000 (/i-iMpc)3. We 



to fc„„x = 0.31 /iMpc ^ ^ with a bin size given by Afc = 
0.015 /iMpc~^ We consider Nk = 20 fc-bins for the 
power spectrum, corresponding to Nt = 1015 triangle 
bins, including all triangle shapes and orientations corre- 
sponding to 7.5 X 10^° elementary triangles . We use 6000 
realizations of the survey |59l l6(]l | , such a large number 
is needed in order to estimate covariance matrices larger 
than 10^ X 10^ in size (see next section). 

For the mock catalogs of the LRG sample, we use 
6000 mock catalogs constructed with PTHalos, using the 
following cosmological parameters: dark matter density 
Q,d — 0.229, baryon density Vti, = 0.046, cosmologi- 
cal constant with density JIa = 0.725, Hubble constant 
h = 0.71 and fluctuation amplitude a% = 0.75 at the 
mean redshift of the survey of Zmcan — 0.35. In these 
mock catalogs the LRG galaxies populate dark matter 
halos according to an HOD prescription |33| for the mean 
number of galaxies in a halo of mass m 



(NsAm)) 



-fl'l'mi-n./l'l't- 



1 



/ ?7l \ ' 

Vmi/ 



(7) 



where the first contribution is that due to a central galaxy 
(with nearest integer scatter), the rest being satellite 
galaxies which are taken with a Poisson distributed scat- 
ter [33. The parameters are chosen by a best fit pro- 
cedure of the larg e-scale redshift-space correlation func- 
tion given in |64| (including the survey covariance ma- 
trix) and the small-scale redshift-space correlation func- 
tion given in p3q |. The resulting parameters are TTimin = 
5x1O13M0/i-\toi = IO^^Mq/i-i anda = 1.95. Given 
Eq. lO, the large-scale bias parameters are given by, 

bi = — I dm n{m) (iVgai(m)) bi{m), (8) 

ng J 

where n(m) is the halo mass function (assumed to be 
that in l36|). bi{m) are the corresponding halo bias pa- 
rameters p37| and the galaxy number density is given by 



dm n{vi) {Ngax{m)). 



(9) 



measure the power spectrum from k„ 



0.02 /iMpc" 



For the parameters given above, n^ = 8 x 10~^, 61 = 2.11, 
62 = 1-1 and 63 = —2.8. In practice, the values of the bias 
parameters measured in the mock catalogs are slightly 
different from the analytical calculations due to the par- 
ticular prescription adopted to describe how individual 
halos are populated. We find, for example, 61 — 2.17 and 
we use this as our fiducial value for the LRG large-scale 
linear bias. The mock catalogs have the radial selection 
function expected by the end of the SDSS survey and de- 
scribed in W^ and an angular selection function equal to 
unity inside the survey region. The redshift-space density 
field in them is weighted according to the FKP method 
with Po = 40, 000 ( /i-i Mpc)3. 

In section IIVI below we discuss in more detail how we 
take into account redshift distortions. In brief, we as- 
sume the 2LPT or PTHalos simulations give the correct 



answer (note that these do not assume perturbation the- 
ory for the real-to-redshift space mapping), and compute 
deviations from the fiducial model by tree-level perturba- 
tion theory. A full discussion about accurate theoretical 
predictions for statistics of galaxies in redshift-space and 
their possible systematics is beyond the scope of this pa- 
per, and will be presented elsewhere. 



III. COVARIANCE MATRICES 

In order to perform a joint likelihood analysis of the 
power spectrum and bispcctrum, as detailed in section llVl 
below, we need to compute their covariance properties. 
The full covariance matrix Cij obtained from our mocks 
catalogs by measuring the power spectrum and bispcc- 
trum, is defined as 



where 



Cij = {6Xi5Xj 



(10) 



where 5Xi = Xi — Xi and Xi equals the power spectrum 
Pi for i — I, ...jNk with Nk the number of power spec- 
trum bins, or the bispectrum for i = N^ + l, ■■■,Nk + Nt, 
with Nt the number of bins in triangle space. 

In what follows, we consider the three contributions 
to the general covariance matrix Cij, that is {SPi6Pj), 
(SBiSBj) and (SPiSBj) in turn, and compare the ex- 
pected contributions to the values measured from the 
mock catalogues. 



A. Power Spectrum Covariance 

Our power spectrum estimator can be written as 



Pik) 



^^d3,,^d3,,^,(q,,).,,^,„ (11) 



TABLE I; Power spectrum cross-correlation coefficients Clj 
between different scales as measured from the main sample 
SDSS mock catalogues. The value of the wavenumber k for 
the corresponding bin is given in the first column in units of 
hMpc~^ . For brevity, only even bins are displayed. 
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Vp(fc) = /d^gi /d^ga SD{qi2) ^ ^Trk^Ak, (12) 

Jk Jk 

and the integral over the bin k of width Afc is given by 



(./c+A/c/2 

d^q= / dfcfc2 

k Jk-Ak/2 



an. 



(13) 



The bin width Ak does not necessarily coincide with the 
fundamental frequency kf (in our analysis we will con- 
sider the case Afc = 3fc/). If the value of the power spec- 
trum averaged over all the realizations is P{k) — (P(fc)), 
it is easy to see that the covariance between power spec- 
trum bins can be expressed as |33, |33 



Cl 



{5P{h)5P{k,)) = 
2 k) , 



Vpih) 



fc3 



+1 



dcos6' T{ki,kj,( 



(14) 



where the first diagonal term is the Gaussian contribu- 
tion and in the second, non-Gaussian, term T(ki, kj, 9) = 
r(ki, — ki, kj, — kj) is the trispectrum of the dark matter 
field and 9 is the angle between the vectors k^ and kj . 

Note that expressions such as Eq. H14() , or the analogue 
ones for the bispectrum and mixed covariance discussed 
in the next sections, assume that the survey window is 
effectively a delta function in Fourier space, and they are 
therefore exact in the case of a periodic box. However, 
they provide a simple estimate of the covariance for more 
generic survey geometries except for the case of the mixed 
power spectrum - bispectrum covariance, (SPiSBj) as we 
will see in section llll (JI below. Of course, in analyzing our 
mock catalogs we do not make any such approximation as 
the survey geometry is properly taken into account by the 
FKP method. The estimators for power spectrum and 
bispectrum will therefore include convolutions with the 
window function and shot noise terms and the covariance 



TABLE II: Same as Table El but for the LRG sample SDSS 
mock catalogues. For brevity, only one bin every three is 
shown. 
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FIG. 2: 



Power spectrum cross-correlation coefficients r^j between different scales for the main (left) and LRG sample (right) 



as measured from the SDSS mock catalogs. Black indicates maximum cross-correlation (r^j = 1) 



white no cross-correlation 



{r(j = 0). The wavenumbers are given in units of Afc ~ 0.015 ft Mpc ^ for the main sample power spectrum, while Afc 
0.0075 /iMpc~^ for the LRG power spectrum. Note that bin 27 in the LRG case corresponds to bin 13 in the main sample 
case: the lower left region in the main sample plot encloses the scales corresponding to the LRG plot. 



will then be computed from the measured statistics in 
each mock catalogue. 

Figure [21 shows the redshift-space power spectrum 
cross-correlation coefficients 



p _ 



ct 



riPnP 



(15) 



for the main (left) and LRG (right) sample power spectra 
measured from our mock catalogs. The values are rang- 
ing from (white) to 1 (black). We used iV^ = 20 bins 
for the main sample, N^^^ — 27 bins in the LRG sample 
case. The numerical value of the cross-correlation coef- 
ficients for the main sample power spectrum is given in 
Table m whereas Table HTl presents the LRG power spec- 
trum case. Note that in this case the maximum value 
for the wavenumber considered is k^^'^ = 0.2/iMpc~ , 
instead of k^£ = 0.3/iMpc~^ for the main sample. 

As evident from these tables and Fig. 12 the cross- 
correlation between different scales is stronger for the 
LRG power spectrum case. For example let's consider for 
the main sample the coefficient C^g = 0.34 (Tab.Hl) cor- 
responding to the wavenumbers 0.126 and 0.63 /iMpc^^ 
and compare it to the LRG coefficient C|^i7 — 0.66 
(Tab. ini corresponding to the wavenumbers 0.123 and 
0.60ft.Mpc^^. If LRG galaxies were simply a linearly 
biased population compared to the main galaxy sample 
(here assumed to be unbiased), then one would expect 



exactly the opposite given our choice of bin widths, that 
is a smaller value in the LRG case. This is so because the 
cross-correlation coefficient is independent of the volume 
of the sample (which appears in Eq. H14|l only through 
kf), proportional to the amount of non-Gaussianity (here 
given by the averaged trispectrum divided by the power 
spectrum squared), and proportional to the bin width 
Afc. The reason for this last dependence is that the non- 
Gaussian noise term does not get beaten down by bin 
averaging, whereas the Gaussian one (which dominates 
in the denominator in Eq. H15|) ) does jSg. Since our bin 
size for LRG power is half that of the main sample, and 
linear bias does not alter non-Gaussianity (and redshift- 
distortions do so very slightly at large scales, and fur- 
thermore cTg is somewhat larger in the main sample) , this 
would give that LRG power should be cross-corrrelated 
about half as much as main sample galaxy power. 

The difference in behavior is thus a reflection of non- 
linearity in the LRG bias, which creates additional 
non-Gaussianity. In fact, this is expected in standard 
scenarios of galaxies, since Eq. ^ naturally predicts 
that for galaxies that populate high-mass halos where 
6i ~ 2, 62,63 should be at least of order unity. How- 
ever, we caution t hat, unlike the case of linear bias 
61 (m) [43, m El, H m, the expressions for non- 
linear bias parameters for halos 62 (w), 63(771) given by the 
peak-background split |37l 143 . 144. |45| (also assumed by 
PTHalos) have not been tested against current numerical 
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FIG. 3: Bispectrum cross- correlation coefficients r^j among the first 19 triangles at large scales (left), and among two sets of 
the 20 triangles at the smallest scales we consider (right). The triplets indicate the wavenumbers of the triangles sides in units 
of the ^-bin width, Ak c^ 0.015 /zMpc"^ 



simulations (see |4y| for early work). This is an impor- 
tant issue since the prediction is that 62,^3 ^^^ strong 
functions of halo mass for the range relevant for LRG 
galaxies (see e.g. Fig. 8 in [33)? and small changes in the 
HOD parameters that leave the linear bias within obser- 
vational bounds can change the non-linear bias param- 
eters significantly. It is for this reason that we do not 
consider the bispectrum of LRG galaxies in this work, 



TABLE III: Bispectrum cross-correlation coefficients for tri- 
angles at the largest scales. Each triangular configuration is 
given in terms of the three vectors ki, k2, ks in units of the 
k-hin width, Ak cri 0.015 /iMpc"^ 
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since its prediction has significant uncertainties. We are 
currently working on addressing these issues. 

The nonlinearity in the bias relation for LRG galaxies 
implies that the power spectrum can only be reasonably 
well approximated by linear bias up to larger scales than 
in the main sample case. This is why we take fc,,,-,^ — 
0.2/iMpc"^ for LRG's, instead of A:„„, = 0.3/iMpc"^ 
for the main sample. These are reasonable, though some- 
what arbitrary values. In practice the allowed values of 
A:„ax can be empirically tested by looking at higher-order 
correlations and looking for scale-dependence in the de- 
rived bias parameters p!M[l9J|. 



B. Bispectrum Covariance 

In an analogous way to the power spectrum case, we 
can define the estimator for the bispectrum, |5lj| . 

B{kuk2,h) = -f fd^qi (d\2 fd'^qs 
Vb Jki Jk2 Jh 

x5D(qi23)'5qi(5q2^q3, (16) 

with 

Vb = jd^qi d\2 d^qsSoiflus) 

~ Stt^ fcifcafca Afc^ (17) 

then the covariance between two triangle configurations 
(where i and j represents the triangles while (ii,i2,*3) 
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FIG. 4: Mixed cross-correlation coeHicients 



between main sample power spectrum and bispectrum. We show the 20 



largest scale triangles (left) and the 21 smallest scale triangles (right) against all power spectrum bins. All numbers indicate 
wavenumbers in units of the fc-bin width, Afc ~ 0.015 hMpc~^ . The vertical bands are due to beat-coupling. 



and (ji,i2,j3) are the corresponding fc- vectors triplets) 

is, 



Cf^ = {6B,SB, 



l3 



VByi') 

x5D(qi23)'5n(qi+P23)S(qi, P2, P3)S(qi, q2, q3) 
+cyc. 



~T~"iljl ' 



k} 



d\i... d\Jd^P2 d^P3 



'''VBit)VBij)jk: 

X^£)(qi23)'5u(p23-qi)^('?l)r(q2,q3,P2,P3) 

+cyc. 

+17TTw^ A'^i . . . fdhs fd^Pi • . . /d^P3 

VB{i)VB{j) Jk,^ Jk,^ Jk,, Jk,^ 

X(5D(qi23)fo(pi23)'?6(qi, q2, q3, Pi, P2, P3), (18) 

with T6(ki, ..., kg) representing the 6-point connected 
correlation function in Fourier space. At large scales, 
the main contribution to the variance of the bispectrum 
is Gaussian and therefore 

^3 

AB^h, k2,h) ~ SBjf- P{h)P{k2)P{h), (19) 

Vb 

with SB = 6, 2, 1 for equilateral, isosceles and general 
triangles, respectively. 



From the expression for Cf, we see that the largest non- 
Gaussian contribution to the extra-diagonal elements of 
the bispectrum covariance matrix should arise in triangu- 
lar configurations sharing two sides, with an extra factor 
when these are equal sides of isosceles triangles. Such 
large terms can be easily identified in Fig. |31 where we 
show the bispectrum cross-correlation coefficients. The 
value of the cross-correlation coefficients at the largest 
scales is given in Table IIIII Note that, even at small 
scales, the bispectrum cross-correlation coefficients re- 
main small, with values usually lower than 0.3, typically 
quite smaller than in the power spectrum case. 



C. Mixed Covariance: Beat Coupling 

Given the estimators for the power spectrum and bis- 
pectrum defined in Eqs. (|ll|l and (|16() . the mixed terms 
in the general covariance matrix are 



Cr,^ ^ (SP.SB,) = 



•J I] I 



2k) 

k) 



P(h)B{kj^,kj^,kj^) + eye. 



'd\i d'p,... d'ps 



Vpil)VBij) Jk 
X^D(pi23)7'5(qi, -qi, Pi, P2, P3 



(20) 



where T5(ki, . . . , ks) stands for the 5-point connected 
correlation function in Fourier space. At large scales. 



the first term in Eq. (|20|l dominates, and moreover, this 
is expected to be an important contribution. To see this, 
compare its magnitude to the expected signal 






2sBk) 



sb_ 

2ti 



k> 



(21) 



which is comparable to the same ratio for the diagonal 
covariance of the power spectrum. 



P? 



2k) 

Vp{i) 



1 

2^ 



(22) 



Figure 0] shows the cross-correlation coefficients between 
the first 20 (left) and last 21 (right) bispectrum config- 
urations and all power spectrum bins. The terms just 
discussed correspond to the diagonal in the right panel, 
and a few of the elements in the bottom part of the left 
panel, where the power is calculated at one of the sides 
of the triangle. The value of the mixed cross-correlation 
coefficients at the largest scales is given in Table Hvl 

However, it is evident that there are significant cor- 
relations beyond these, for triangles which include the 
smallest value of /c as a side with every bin of the power 
spectrum. Indeed, Eq. (|20|l ignores important contribu- 
tions that dominate the mixed covariance matrix. The 
reason is that so far we have ignored the effects of the 
window of the survey. 

In a finite survey of size ~ L, the uncertainty principle 
implies one cannot measure Fourier modes to a better 
accuracy than 5k ~ tt/L, since two waves of frequency k 
and k±Sk differ only by half an oscillation from one end 
to the other of the survey, i. e. there is not enough room 
inside the survey to tell them apart. This implies that in 
reality the power spectrum estimator in Eq. (|11|) written 
in terms of the observed Fourier modes will necessarily 
contain, due to the survey window, cross-terms in the 
underlying Fourier modes, written schematically as 



^q ^-q+e, 



(23) 



TABLE IV: Mixed cross-correlation coefficients between 
SDSS main sample bispectrum and power spectrum at the 
largest scales. Each triangular configuration is given in terms 
of the three vectors fci, ^2, fca in units of the fc-bin width, 
Afc~ 0.015 ft Mpc"^ 
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where e < Sk, apart from "true power" contributions 
(5q(5_q — |(5qp. Although such terms do not contribute to 
the expectation value of the power, they do correlate very 
well with appropriate bispectrum configurations. Indeed, 
due to quadratic nonlinearities two nearby Fourier modes 
couple to the beat mode between them. 



^q <^-q+e ^ <5q ^2(-q, e) S-q Se 



(24) 



which means that these terms dominate the fluctuation 
in power at high wavenumbers where g 3> e, giving the 
non-intuitive result that the errors of the power spectrum 
in the nonlinear regime are dominated by the large-scale 
power j^ 1^ . From Eq. (|24|l it follows that such terms 
cross-correlate very well with the bispectrum of isosceles 
triangles with one small side of the order of the survey 
window ~ e. 



{Sk S-k^ 



Sp 6^p (5e 



P{k)Pip) P{e). 



(25) 



Therefore, for all k we expect power spectra to cross- 
correlate with bispectra of "narrow" isosceles triangles. 
These are the vertical features seen in Fig. 0] 

Beat coupling implies that the whole power spectrum 
and the bispectrum of "narrow" isosceles triangles fluctu- 
ate together depending on the large-scale power. As we 
shall see in section rV(]l this has interesting implications 
for the likelihood analysis. 



IV. THE LIKELIHOOD FUNCTIONS 

We now consider a hypothetical joint analysis of large- 
scale structure (LSS) and cosmic microwave background 
(CMB) anisotropics. In order to be specific and illustrate 
the amount of information that we expect to extract in 
the very near future, we consider the first year WMAP 
data, the power spectrum and bispectrum of the SDSS 
main sample of galaxies, and also include the SDSS power 
spectrum of the luminous red galaxies (LRG). The SDSS 
"data" is obtained from the mock catalogs described in 
section HTl and corresponds to the survey in its expected 
final form. In this section we describe the LSS and CMB 
likelihood functions that we use to derived the constraints 
discussed in the next sections. 



A. The LSS likelihood 

For simplicity we assume that the power spectrum and 
bispectrum estimators are Gaussian distributed. This 
is certainly a good approximation near the maximum 
wavenumber we include, but becomes worse at large 
scales, where only a few modes (for the power spec- 
trum) or triangles (for the bispectrum) contribute. The 
deviations from Gaussian likelihood can be included by 
calculating the likelihood from the Monte Carlo pool of 
mock catalogs 53]. Ignoring the non-Gaussianity of the 



likelihood can lead to a biased estimation of parame- 
ters [53, I53 • Since here we are only trying to understand 
the improvement brought by adding the bispectrum to 
the standard joint analysis of CMB and LSS, and most 
of the information added by the bispectrum is coming 
from scales small compared to the survey, our assump- 
tion should be safe. The combined power spectrum and 
bispectrum likelihood function is then 



where 



In £ = In Cp + In Cb + In CpB , 






Nt Nt 

i=i 3=1 



and 



ln£ps = 



j=i 3=1 



SB, 



(26) 



(27) 



(28) 



(29) 



takes into account the mixed elements of the inverse 
covariance matrix C^^. In Eqs. (|27I29|) . the indices i 
and J run over the bins in /c-space for the power spec- 
tra, Nk in all, as well as over the Nt configurations in- 
cluded in the bispectrum analysis. Also, 5P = P, — P* 
and SB = Bs — B*, where Ps = Ps(p;fc) and Bs = 
Bs{p; /ci, A:2, /ca) are the redshift space galaxy power spec- 
trum and bispectrum as a function of the parameters p 
while P* = Ps{p*;k) and B* = Bs{p*;ki,k2, fcs) are the 
redshift space galaxy power spectrum and bispectrum of 
the fiducial model (with parameters p*). 
In the most general case we consider 



{T,As,uJd,uJb,^A,ns,w), 



(30) 



defined as the reionization optical depth, r, the primor- 
dial amplitude of scalar fluctuations, As, the physical 
dark matter density, ujd = ^dh^ , the physical baryon 
density, uji, = 0.bh? , the dark energy density, 0.^, the 
scalar spectral index, Us and the dark energy equation 
of state parameter, w = p^l Ph- The bias parameters in- 
clude the main sample linear and quadratic bias, 5i and 
62, and the LRG linear bias h^^'^ . 

The covariance matrices are calculated at maximum 
likelihood from our mock catalogs, that is, we do not 
include a possible dependence on parameters to be esti- 
mated. A simple, approximate, check of such dependence 
on the bias parameters did not yield appreciable differ- 
ences in the final results we present in section Ivl 

When we study below results from the power spec- 
trum or bispectrum individually the inverse matrix C~^ 
in Eq. (|27|l or H28|l will be replaced by the inverse of the 
individual matrix C/J 



will study as well the case of combining the two statis- 
tics without taking into account their mixed covariance, 
in which case also only C^ and C^ will be needed. 

The likelihood function for the LRG power spectrum 
is equivalent to Eq. H27|l and the corresponding covari- 
ance matrix is independently determined from the LRG 
mock catalogues. Since the mean redshift of the LRG 
sample is 2: ~ 0.35 compared to 2; ~ 0.1 for the main 
sample, with little overlap, we assume the two samples 
are independent. 



1. Power Spectrum 

Deviations from the fiducial redshift space power spec- 
trum monopole P*, as a function of the parameters p, 
are modeled in the following way 



P«(p;fc)=6?/f(p)/(p;A:)P;(fc), 



(31) 



where P* is the power spectrum measured from our 
mocks catalogs in redshift space and where 



/(p;fc) 



A 
]4* 



r(p;fc) 



T*{k) 



k_ 
k^ 



Dip)/n, 



(/i*)";' 



(32) 



r(p; fc) being the transfer function, D+{p) the growth 
factor and kp = 0.05 Mpc~^ the pivot point correspond- 
ing to the scale whose power is unaffected by varying the 
spectral index, and 



fP/ ^ _ ao(/3) 



/r(p) ^ 



aoiP* 



(33) 



(34) 



{SP.SPj) or Cf. = {SB,SBj). We 



is the redshift-space correction where 

2 1 , 

ao = 1 + 3/3 + -/3^ 

with /9 ~ QJ. /6i, corresponds to the power spectrum 
monopole [5g. Note that we use Eq. (|34|l only to model 
deviations from our fiducial cosmology assumed in the 
mock catalogs that include nonlinear effects from the 
redshift-space mapping. We assume a fiducial model that 
is unbiased, 61 = 1 and 62 — 0. The likelihood function 
for the LRG power spectrum is computed in the same 
way, except for the fiducial value of the LRG linear bias 
parameter h{^'~' = 2.17. 

Note that, since redshift distorsions break the statis- 
tical isotropy expected in real-space, the redshift-space 
power spectrum is a function of the direction as well as 
the magnitude of the wavevector k. In this work we in- 
clude only the monopole of the redshift-space power spec- 
trum, i.e. the average over the angle formed by k and the 
line of sight. In principle one can take advantage of the 
angle-dependence by measuring the quadrupole term in 
the Legendre polynomial expansion and obtain a better 



10 



determination of the /3 parameter, further strengthening 
the constraints presented in section IVI below. 

We calculate the transfer functions from CMB- 
FAST ^0], computing the value of T(p; k) for every value 
on a limited grid in parameter space, then interpolat- 
ing over the final parameter grid for each value of the 
wavenumber k involved in the analysis. 

For the growth factor £'(17™, w) we take advantage of 
the fitting formula provided in 61]. This is given by 



D{a)~aexp< / din 



a[n^{ay^l]\, (35) 



with 7 — 0.55 + 0.05(1 + w) and where a is the cosmo- 
logical scale factor. 



2. Bispectrum 

We describe the deviations of the bispectrum from the 
fiducial model of the mock catalogs by means of Eulerian 
perturbation theory. Since we are averaging the redshift- 
space bispectrum over triangles with all possible orien- 
tations, similarly to the power spectrum case, we only 
need the monopole term in a Legendre expansion. We 
will consider the following approximation, see Eqs. (20- 
28) in 13 

B(o) ~ a^(/3)|-F2(ki,k2)P3(fci)P<,(fc2) 



6i 



+af(/3)-fPg(fci)Pg(fc2)+cyc., 



(36) 



where Pg{k) ~ h\P(k) is the galaxy power spectrum, and 



a^(/3)^ 1 + ^/3 +i/?2, 



(37) 



describes the bispectrum monopole redshift space correc- 
tion, obtained from Eq. (24) and (28) in jsj by averaging 
over the angle between ki and k2 and dropping the de- 
pendence on the second-order velocity kernel and velocity 
dispersion which should partially cancel at large scales, 
to approximate the configuration dependence found in 
simulations for the redshift-space bispectrum in |53 |. 

To compute the dependence of the bispectrum on cos- 
mological parameters we will therefore use 

Bsiv) = /f(p,)6?s|j 

+/f (p) 62 hi S, (38) 

where B* is the redshift-space bispectrum measured from 
the mock catalogs and 

B = 2F2(ki,k2)F;(fci)/iF;(fc2)/2+cyc., (39) 
S - P;(fci)./iP;(fc2)/2 + cyc., (40) 



while 



/f(p,P*) 



_ aim 



(41) 



with (3* as fiducial (3 and /i = /(p; fci) as defined above. 



3. Inverting the Covariance Matrix 

The values of the entries of the complete covariance 
matrix dj = (XiXj) with Xi — Pi,Bi span several 
orders of magnitude and thus a direct computation of 
its inverse is susceptible to numerical instabilities. We 
therefore "normalize" the covariance matrix by factor- 
ing out in the X = P, B vector the power spectrum and 
bispectrum predicted by linear theory and Eulerian PT, 
respectively. The resultant entries for this "normalized" 
covariance matrix are therefore all of order unity. Still, 
by performing a singular value decomposition (SVD) one 
can notice a poor determination of a few singular values, 
about 17 out of 1035 in the complete, k^i^ = 0.3 h Mpc~^ 
case, indicating that 6000 mock catalogs is enough to de- 
termine most of the elements except for a small fraction. 
In the final analysis we compute the inverse by means of 
its SVD inverse by dropping these few singular values, 
assuming this might be a sign of a not optimal deter- 
mination of the matrix Cij . By doing this we make a 
conservative choice since the operation amounts to dis- 
card some of the potential information contained in the 
covariance matrix. Therefore our final error bars increase 
slightly. In computing the inverse of the individual power 
spectrum and bispectrum covariance matrices no such 
limitation is needed. 



B. The CMB likelihood 

To combine the results of the LSS likelihood analy- 
sis with CMB data as measured by WMAP 1-year, we 
need to compute the CMB anisotropics power spectrum 
and its corresponding likelihood for each model in our 
grid. This procedure is computationally expensive: ver- 
sion 4.5.1 of the CMBFAST code 0| takes about 30 
seconds per model, and the WMAP likelihood code for 
the first year data release [43 takes 2 — 4 seconds. 

A possible a ppr oach to reduce the computing time is 
investigated in |57| , where a polynomial approximation to 
the multidimensional log-likelihood is computed, allow- 
ing for an evaluation of the likelihood of each model in 
tenths of a second. However, their available code, CMB- 
FIT 1.0, does not include the dark energy equation of 
state parameter, w. Although calculating such a poly- 
nomial fit still requires sampling the likelihood surface, 
it has to be done only once, thus reducing enormously 
the computational time needed for any ulterior likelihood 
analysis. 

Motivated by that idea, we compute a polynomial fit 
to the CMB likelihood function based on the 7 param- 
eters in Eq. H30|l . We use uniform priors in the follow- 
ing ranges: < r < 0.3, 0.018 < cob < 0.028,-2 < 
w < 0.2,0.5 < ^ < 1.4,0.4 < r^A < 0.9,0.08 < ujd < 
0.22,0.8 < Us < 1.1. 

We compute the likelihood on a homogeneous grid 
with 15 points per dimension for the parameters 
(wrf, Wh, 17A,n3, w), and 30 points for {t,As). To speed 
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up the calculation, we divided the problem in two steps: 
firstly, we computed the likelihood on the grid approxi- 
mating the dependence of the power spectrum on r with 
the multiplicative factor e^^'^. Out of these 7 x 10* ap- 
proximate likelihood values we selected a connected sub- 
set of 263,022 models containing the maximum and de- 
fined by a threshold chosen to be 10 orders of magnitude 
smaller than the maximum. We then recomputed the 
likelihood for the reduced subset with the correct t de- 
pendence. 

We fitted a 4th order polynomial to the log-likelihood 
surface spanned by our reduced dataset using a weighted 
least squares method. We weighted the fitting error of 
each model with its likelihood to counterbalance the fact 
that our grid was relatively coarse and there were many 
more low likelihood models than high likelihood ones. 
The covariance matrix of our 7-dimensional reduced set 
of models is given by C^ = {piPj) - {pi){pj), with (pi) = 
"^PiClp). In order to improve the numerical behavior 
of the fitting algorithm we first changed from p— space 
to the variables z with zero mean and unit covariance 
defined as 



z = E(p-(p)), 



(42) 



where the rows of E were defined as the eigenvectors of 
the covariance matrix C divided by the square root of 
their corresponding eigenvalues, i.e. such that ECE* — 
I, and thus (zz*) — I, 57). The 4th order polynomial, 
containing M — 330 terms, can be written in terms of 
the new variables as 



y = \nC = qo 



7 

41 = 1 



?1 
9l 



7 

12 — ii 



r 'i-i'i- 
[l2 



7 

E 



^^3 Us 



7 

E 

44=13 



Qi 



(43) 



In order to make this expression compact we arranged 
all possible products of z^'s (up to the 4th power) into an 
M-dimensional array x. 



X = {1, 21, ... , Z7, ZlZl,ZlZ2, . . . , Zj}, 

and the corresponding coefficients q into, 

n - Sn «1 n'^ n^^ «7777l 

so that Eq. l|I5|) could be casted as 
y = X • q. 



(44) 
(45) 
(46) 



Next, we arranged our N = 263,022 datapoints, written 
in the x format, into an {N x M) matrix X and their 
corresponding likelihoods into an A'^-dimensional vector 
y. Therefore, the weighted error of the polynomial fit 
was given by. 



1 ^'' 
mE"'^(2^» 



M 

E 

j=i 



^ijlj) 



1/2 



(47) 




FIG. 5; WMAP TT-I-TE marginalized likelihoods for the 7 
cosmological parameters in our wCDM model. Dots corre- 
spond to the marginalization over the reduced dataset in the 
grid (263,022 points) while the solid line was obtained using 
the weighted polynomial fit. 



where Wi = C{pi). The coefficients qi were then chosen 
such as to minimize e, 



q=(X*X: 



r'x'y, 



(48) 



where we defined X, 



w.Xi 



In order to avoid unphysical likelihood values due to 
polynomial artifacts in low confidence regions that were 
poorly sampled, we replaced the polynomial fit by a sim- 
ple Gaussian for z — |z| outside the 2-sigma level, given 
that the likelihood distribution has a spherically symmet- 
ric tail for z > 2.5, well approximated by £ ex e~^ '^. 

To test our weighted fit we compared it, for w = 
-1, against CMBFIT for a 6-parameter ACDM model 
(t, As, tt'd,a;f,, Oa,'^s), finding very good agreement. An- 
other estimator of the goodness of the fit was the rms 
fitting error (A In/:). We found that (Aln£) = 0.29 for 
our fitting dataset of 263,022 models. Furthermore, we 
made a Monte Carlo Markov Chain test of 3000 models, 
which yielded (Aln£) = 0.32. These errors are similar 
to those reported in |57| for the 7-parameter models of 
CMBFIT. 

Figure shows the first year WMAP TT-hTE 
marginalized likelihoods for the case of the 7-parameter 
wCDM models, Eq. l|Sni) . obtained by a grid marginal- 
ization over the reduced dataset (dots) versus the ones 
using the weighted polynomial fit (solid line) . This shows 
the procedure described above is robust enough for our 
purposes. 
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TABLE V: Fiducial values for the cosmological and bias pa- 
rameters assumed. 



Parameter 


Fiducial value 


(^d 


physical dark matter density 


0.1222 


UJt 


physical baryon density 


0.0232 


Oa 


dark energy density 


0.699 


Us 


scalar spectral index 


0.977 


As 


scalar fluctuation amplitude 


0.81 


w 


dark energy equation of state 


-1 


T 


reionization optical depth 


0.124 


bi 


main sample linear galaxy bias 


1 


b2 


main sample quadratic galaxy bias 





f^LRG 


LRG linear galaxy bias 


2.17 


Derived parameters 


Fiducial value 


0-8 


galaxy fluctuation amplitude 


0.917 


fim 


matter density 


0.301 


fit 


baryon density 


0.048 


h 


Hubble parameter 


0.695 



RESULTS 



In this section we present the results of the hkelihood 
analysis in two classes of fiat cosmological models. The 
first, section l\r^ corresponds to ACDM models depend- 
ing on six cosmological plus three bias parameters: the 
density parameters lu^, oJb, ^^a, the spectral index ng, the 
fluctuations amplitude As , the reionization optical depth 
T plus the linear and quadratic galaxy bias coefficients 
6i and 62 for the main sample and the linear bias for 
the LRG sample, hi^^ . In the second class, denoted as 
wCDM models, section IV El we allow for a dark energy 
equation of state parametrized by the ratio of pressure 
to energy density w, assumed to be constant. 

We include the temperature and polarization WMAP 
1-year likelihood by means of the interpolation fit de- 
scribed in section llVBI ffor an update to the 3-year data, 
see the Appendix) . Wc introduce here a flat prior on r by 
limiting its values from zero to 0.3. The difference with 
the case of t taking values up to 0.8 is negligible (tested 
for w = —1) and, most importantly, such a prior is more 
than justified by the three-year WMAP data which favors 
values of r close to 0.1 |3.l5a|. 

The fiducial values chosen for the present analysis are 
given in Table IVl Note that they do not coincide with the 
maximum likelihood values obtained from the WMAP 
data alone, rather they correspond to those obtained for 
the WMAP-I-SDSS power spectrum 6 parameters case 
in [J. These values are only relevant in the sense that 
they determine the point in parameter space about which 
we compute the errors. As long as this point is realistic, 
the results we present should be insensitive to their pre- 
cise values. 



A. ACDM models 

We present now the expected errors on cosmological 
parameters from an analysis that considers different com- 
binations of the main sample power spectrum and bis- 
pectrum and the power spectrum of the LRG sample 
with WMAP CMB data. We restrict here to the case 
of ACDM models, i.e. w — —\. The results for the 1-cr 
marginalized uncertainties are given in Table IVII where 
we show, for clarity, the average between upper and lower 
limits. 

To see more clearly the benefits brought by using dif- 
ferent statistics, in parenthesis we indicate the fractional 
improvement over the WMAP plus main sample power 
spectrum case (W-l-P), defined as 



improvement factor 



iW+P 

A 



1, 



(49) 



so a 50% (100%) improvement corresponds to reducing 
the errors by a factor of 1.5 (2). 

The first two columns in Table IVll show that analyzing 
the power spectrum and bispectrum separately can pro- 
vide similar constraints (with the bispectrum determin- 
ing an extra parameter, 62). This can provide important 
consistency checks, as the Gaussian and non-Gaussian 
properties of galaxy clustering must yield consistent re- 
sults. 

As expected, the effectiveness of the bispectrum in con- 
straining cosmology depends significantly on the small- 
est scale considered due to the fast rise in the num- 
ber of triangles available. One can notice how, already 
at fc^^f — 0.2/iMpc~ when combined with the power 
spectrum, it can improve errors by a 13 to 56%. At 
k^^l = 0.3/iMpc~^ when considered alone with CMB 
information the bispectrum can actually improve over 
the power spectrum by 6% to more than 70% for f^A, 
although at the expense of a poorer determination of the 
linear bias. We should keep in mind that the bispectrum 
analysis introduces an extra parameter, the quadratic 
bias, and that one can expect a better constrain on the 
linear bias when combined with the power spectrum. 

A quick glance at Table IVll shows that most of the im- 
provement (numbers in bold) brought by the bispectrum 
are in parameters related to the overall amplitude of fluc- 
tuations and the effective spectral index. This is expected 
as the bispectrum breaks the degeneracy between bias 
and dark matter amplitude fluctuations [a, |l3| , and its 
configuration dependence is sensitive to the spectral in- 
dex because of the anisotropy of tidal gravitational fields 
and velocity flows [63 ■ 

In Fig. we compare the CMB power spectrum likeli- 
hoods to the combined power spectrum, bispectrum and 
LRG power spectrum likelihoods. From this and Ta- 
ble IVII one sees that most of improvement over CMB 
alone is coming from pairs of statistics that involve the 
bispectrum {i.e. either P+B or P^-l-B, not shown for 
clarity) . This is so because the most significant improve- 
ments arise due to breaking of degeneracies present in 
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TABLE VI: ACDM models: expected marginalized errors (68% CL) for WMAPl (temperature and polarization, W) combined 
with the SDSS main sample power spectrum (P) and bispectrum (B) and with the LRG power spectrum (Pl). The percentage 
in parenthesis indicates the improvement over the analysis including the main sample power spectrum alone (W+P), numbers 
in bold indicate errors down by at least 1.5. In brackets we quote the W+P+B errors obtained by ignoring the mixed power 
spectrum - bispectrum covariance. 





W+P 


W+B 


W+P+B 


VV+P+B (no mix. 


cov.) W+P+B+Pl 


W+P+Pl 


W+Pl 




kt'I = 


].2/iMpc-^ 








kt^f =0.2/iMpc-^ 


Auja 


0.0035 


0.0041 (-15%) 


0.0031 (13%) 


0.0030 (17%) ] 


0.0025 (40%) 


0.0026 (35%) 


0.0029 (21%) 


Aujt 


0.00093 


0.00098 (-5%) 


0.00078 (19%) 


0.00082 (13%) ] 


0.00074 (26%) 


0.00081 (15%) 


0.00087 (7%) 


A^A 


0.0133 


0.0113 (18%) 


0.0085 (56%) 


0.0085 (56%) ] 


0.0063 (111%) 


0.0078 (70%) 


0.0091 (46%) 


Aus 


0.022 


0.024 (-8%) 


0.0158 (39%) 


0.0176 (25%) ] 


0.0140 (57%) 


0.0171 (28%) 


0.020 (10%) 


AAs 


0.091 


0.094 (-3%) 


0.064 (42%) 


0.074 (23%) ] 


0.062 (47%) 


0.078 (17%) 


0.085 (7%) 


At 


0.052 


0.053 (-2%) 


0.039 (33%) 


0.044 (18%) ] 


0.038 (37%) 


0.047 (11%) 


0.049 (6%) 


Abi 


0.086 


0.113 (-24%) 


0.060 (43%) 


0.074 (16%) ] 


0.054 (59%) 


0.070 (23%) 


- 


Afea 


- 


0.069 


0.054 


0.062 ] 


0.051 


- 


- 


^^LRG 


- 


- 


- 


0.099 


0.125 


0.137 


Ao-8 


0.068 


0.074 (-8%) 


0.047 (45%) 


0.054 (26%) ] 


0.043 (58%) 


0.054 (26%) 


0.061 (11%) 


Ah 


0.0152 


0.0140 (9%) 


0.0101 (50%) 


0.0108 (41%) ] 


0.0087 (75%) 


0.0106 (43%) 


0.0127 (19%) 


Aftt 


0.00151 


0.00141 (7%) 


0.00124 (22%) 


0.00124 (22%) ] 


0.00111 (36%) 


0.00117 (29%) 


0.00127 (19%) 




fc*lf = 


].3/iMpc-^ 




kif^ ^0.2hMpc-^ 


AcUd 


0.0033 


0.0031 (6%) 


0.0029 (14%) 


0.0026 (27%) ] 


0.0024 (37%) 


0.0026 (27%) 


Acoi 


0.00090 


0.00083 (8%) 


0.00073 (23%) 


0.00072 (25%) ] 


0.00070 (28%) 


0.00080 (12%) 


A^A 


0.0112 


0.0065 (72%) 


0.0063 (78%) 


0.0057 (96%) ] 


0.0052 (115%) 


0.0073 (53%) 


Aus 


0.021 


0.018 (17%) 


0.014 (50%) 


0.013 (61%) ] 


0.012 (75%) 


0.016 (31%) 


AAs 


0.087 


0.081 (7%) 


0.053 (64%) 


0.063 (38%) ] 


0.052 (67%) 


0.077 (13%) 


At 


0.050 


0.047 (6%) 


0.033 (52%) 


0.039 (28%) ] 


0.033 (52%) 


0.046 (9%) 


Afei 


0.081 


0.094 (-14%) 


0.051 (59%) 


0.060 (35%) ] 


0.046 (76%) 


0.068 (19%) 


A62 


- 


0.045 


0.041 


0.044 ] 


0.040 


- 


^l^LRG 


- 


- 


- 


0.084 


0.123 


Ao-8 


0.064 


0.059 (8%) 


0.037 (73%) 


0.041 (56%) ] 


0.034 (88%) 


0.053 (21%) 


Ah 


0.0132 


0.0095 (39%) 


0.0082 (61%) 


0.0080 (65%) ] 


0.0082 (61%) 


0.0101 (31%) 


AQb 


0.00138 


0.00112 (23%) 


0.00111 (24%) 


0.00106 (30%) ] 


0.00104 (33%) 


0.00115 (20%) 



the LSS or CMB ,66]. This manifests itself in the en- 
tries in Table IVll (k'^'^^ = 0.2/iMpc'^) in several ways 
: 1) W+Pl improves mildly over W+P, but W+P+Pl 
improves sig nificantly over W+P (consistent with Ta- 
ble 2 in lill) 2) W+P+B is better than W+P+Pl in 
most parameters (except those related to flm' ^d, ^b 
and JIa, since a better detection of the acoustic scale in 
the LRG sample gives a high quality constraint |64l|). 
This holds even though the signal to noise in B (for 
fc^f = 0.2/iMpc~^) is not as large as in P^ (e.g. com- 
pare W+B vs. W+Pl), because B is more complemen- 
tary than Pi to W+P, i.e. using non-Gaussian infor- 
mation provides a substantially different direction in pa- 
rameter space. When using information up to fc^^f = 
0.3 ft.Mpc""'^, W+P+B constrains all parameters better 
than W+P+Pl except for 104. In this case, adding Pl 
to W+P+B still helps in improving parameters slightly 
(see Fig. IB)), particularly for f^A (or VLm = 1 — ^^a)- 



It is interesting to compare the results on bias param- 
eters to those in a fixed cosmology, as assumed in past 
work [13, [iM H^ HQI • Performing an analysis of the bis- 
pectrum alone with fixed cosmology, one finds for linear 



and quadratic bias the errors (fc, 



MS 



0.3/iMpc"') 



A61 = 0.015, 



A62 = 0.045. 



(50) 



Comparing this to the corresponding entry (W+B) in 
Table IVll we see that when cosmology is allowed to vary, 
the determination of bi suffers from the degeneracy with 
As and, when combined with CMB data, with the optical 
depth T, while the result for 62 is essentially unaffected. 
On the other hand, this is the price one pays for con- 
straining cosmological parameters more accurately. 

Figure 13 shows the marginalized 95% CL contour plots 
of pairs of parameters. The role played by the bispectrum 
in lifting the degeneracy between the galaxy bias param- 
eter 61 and the parameter As determining the amplitude 
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FIG. 6: Marginalized likelihood functions for the ACDM models assuming fc^^x = 0.3/iMpc~^. For the LRG linear bias 
parameter bf only, the dashed line denotes the likelihood obtained from WMAPl plus LRG power spectrum instead of main 
sample power spectrum. 



of dark matter fluctuations is particularly evident. It 
is clear in particular from the bi-Ag contours, that the 
combination of power spectrum and bispectrum, by nar- 
rowing the uncertainty on these two parameters, affects 
the determinations of all the others. The question then 
arises, are the improvements on cosmological parameters 
brought by using the bispectrum just a result of having 
constrained galaxy bias? 



B. Not Just Galaxy Bias 

In order to answer this question, we present in Ta- 
ble IVllI a couple of tests that compare W-l-P with 
W-hP-hB for k^^f = 0.3/iMpc"\ The first two columns 
repeat the constraints shown before in Table IVTl whereas 
the third column shows the W-l-P results when a prior 
on 6i is added to mimic the W-l-P-l-B constraint on bias. 
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FIG. 7: Marginalized 95% contour plots for pairs of cosmological parameters in the ACDM models assuming k^.^^ 
0.3/iMpc~\ 



Since the marginalized likelihood of 6i is approximately 
Gaussian (see Fig. [SJ we can add a Gaussian prior with 
width a given by 



1 1 



'B "P 



(51) 



where as is the error on 6i from the W+P+B analy- 
sis and ap is that from W+P. We see from Table IVIII 
that this reproduces the W+P+B bias constraint closely 



enough. By comparing the rest of the entries in W+P+B 
against W+P+&1 prior it follows that the improvement 
on cosm.ological parameter determination from the bis- 
pectrum is not only due to constraining galaxy bias. 

The right side of Table lVIII presents another test, where 
the bias parameters are fixed (6i = 1, 62 = 0). Compar- 
ing these last two columns we see a significant improve- 
ment from adding bispectrum information. 

The fourth column in Table IVIII shows the analysis 
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TABLE VII: Comparison between W+P and W+P+B with and without galaxy bias assumptions. Includes W+P with a prior 
on 61 to reproduce the error from W+P+B (left) and a case where the bias parameters are fixed (right). The percentage in 



parenthesis indicate the improvement over the respective W+P result. Assumes a ACDM cosmology and k„ 



O.S/iMpc" 





W+P 


W+P+B 


W+P+ fei prior 


W+P+61 prior (Q) 


W+P fixed bias 


W+P+B fixed bias 


Acjd 


0.0033 


0.0029 (14%) 


0.0032 (3%) 


0.0031 (6%) 


0.0031 


0.0026 (19%) 


A(J6 


0.00090 


0.00073 (23%) 


0.00078 (15%) 


0.00073 (23%) 


0.00068 


0.00065 (5%) 


AOa 


0.0112 


0.0063 (78%) 


0.0109 (3%) 


0.0107 (5%) 


0.0106 


0.0062 (71%) 


An^ 


0.021 


0.014 (50%) 


0.016 (31%) 


0.014 (50%) 


0.0112 


0.0096 (17%) 


AA, 


0.087 


0.053 (64%) 


0.061 (43%) 


0.049 (77%) 


0.033 


0.029 (14%) 


At 


0.050 


0.033 (51%) 


0.038 (32%) 


0.032 (56%) 


0.026 


0.024 (8%) 


Afei 


0.081 


0.051 (59%) 


0.051 (59%) 


0.036 (125%) 


- 


- 


A62 


- 


0.041 


- 


- 


- 


- 


Ao-8 


0.064 


0.036 (78%) 


0.043 (49%) 


0.032 (100%) 


0.016 


0.011 (45%) 


Ah 


0.0139 


0.0088 (58%) 


0.0116 (20%) 


0.0110 (26%) 


0.0103 


0.0074 (39%) 


Afit 


0.00139 


0.00111 (25%) 


0.00134 (4%) 


0.00134 (4%) 


0.00134 


0.00104 (29%) 



of the W+P case with a prior on hnear bias, A61 = 
0.036 60], corresponding to the case where the bispec- 
trum is analyzed through the hierarchical amplitude Q 
[see Eq. ©] for a fixed cosmology. We see that in this 
case some of the constraints agree, but the error on bias 
is significantly underestimated, whereas the errors on 17^ 
and h are significantly overestimated. Interestingly, the 
uncertainty in a% is robust to this analysis (which is in- 
correct due to neglecting cross-correlations between Q 
and P and bias with cosmology). 



C. The Effects of Beat Coupling 

It is interesting to see what happens with the con- 
straints on parameters if the mixed covariance between 
power spectrum and bispectrum is ignored, this is given 
in brackets in the fourth column of Table IVII for the 
W+P+B case. Naively, one would expect that excluding 
the mixed covariance should lead to better constraints, 
but as shown in Table IVII this is incorrect for most pa- 
rameters: this is due to the effects of beat coupling. 

As discussed in section llil CI beat coupling means that 
the structure of the mixed covariance matrix is dom- 
inated by up and down correlated fluctuations of the 
whole power spectrum and bispectrum of narrow isosce- 
les triangles depending on the power of the largest mode 
in the survey, as shown in Eq. (|25|) . Not allowing for 
such effect in the covariance matrix means that these 
fluctuations will be mistaken as a signature of larger er- 
rors in the parameters that characterize the amplitude of 
galaxy fluctuations since these are the parameters that 
can mimic such behaviour. Indeed, as seen by compar- 
ing the third and fourth columns in Table IVll including 
the mixed covariance (thus allowing for beat coupling) 
reduces the errors mostly on A^, r, &i, 62 and thus a^. 

In Fig. |S1 we illustrate this point further by showing 
the marginalized 95% contour plots for the linear bias 
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FIG. 8: Marginalized 95% contour plots for the linear bias 
parameter 61 and the quadratic bias parameter 62 in the 
W+P+B case with fe^J = 0.3 /iMpc"^ The top panel shows 
the full analysis with (blue, continuous line) and without (red, 
dashed line) the mixed covariance matrix, whereas in the bot- 
tom panel the two analyses are repeated excluding the lowest 
power spectrum bin, thus suppressing beat coupling. 



parameter h\ and the quadratic bias parameter &2 in the 
W+P+B case with A:;^f = 0.3/iMpc"^ In the top panel 
we show the full analysis, whereas the bottom panel drops 
the bin with the lowest value of k in the power spectrum. 
The top panel shows a significant difference between in- 
cluding the mixed terms (SPiSBj) in the complete co- 



B/Bbbks 




FIG. 9: Baryon acoustic oscillations imprinted in the bis- 
pectrum as a function of fei and ^2 for fixed angle given by 
fci • fe = —1/2. The diagonal fci = fe corresponds to equilat- 
eral triangles. As the angle between the two vectors is varied 
the pattern of peaks moves accordingly. 



variance matrix (solid) and dropping them (dashed) . We 
see that including the mixed covariance gives clearly a 
tighter constraint on the two parameters together with a 
slight degeneracy. 

In the lower panel the same contours are plotted but 
now the analysis excludes the power spectrum bin cor- 
responding to the smallest fc-value, thus suppressing the 
effects of beat coupling. We see that in this case there is 
not much difference between including or not the mixed 
covariance, in fact the mild degeneracy induced by in- 
cluding the mixed covariance leads to slightly larger er- 
rors for 5i and 62. The same behavior is seen with all 
other parameters when the first /c-bin is excluded, ex- 
cept for As and r which still show a minor improvement 
when the mixed covariance is included. This is likely due 
to residual beat coupling, e.g. a careful look at the left 
panel in Fig. ^ shows that vertical features also exist for 
modes with k = 2Afc, although at a much lower ampli- 
tude. 
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square of the power spectrum, the 15% modulation in 
power leads to a 30% modulation in the bispectrum. At 
k — 0.1 hMpc^^ the signal to noise in the bispectrum is 
about twice smaller than for the power spectrum 60], so 
this scale roughly presents the limit after which the bis- 
pectrum gives a better constraint on acoustic oscillations 
than the power spectrum. Unfortunately, at z ~ the 
acoustic oscillations are washed out by nonlinearities for 

k>o.ihMpc~^ liilziliilzi- 

A fair assessment of the improvement on the detection 
of acoustic oscillations by using the bispectrum is beyond 
the scope of this paper and will be presented elsewhere. 
Here we note that our mock catalogs, although not exact 
in their nonlinear properties, do include the suppression 
of acoustic features, that is, we do not assume Eulerian 
second-order perturbation theory as done in Fig. |51 for 
illustrative purposes. 

In order to assess the impact of acoustic features in 
our study we compute marginalized likelihoods using in- 
stead the BBKS fitting formula for the transfer function, 
the results are shown in Fig. 1101 where we reproduce the 
same marginalized 95% CL contour plots given in Fig.Q 
Since this transfer function depends exclusively on the 
shape parameter T = r2m/iexp[— 0^(1 -I- \/2h/flrn)] |75l |. 
we generically expect an enhanced degeneracy between 
rim (or, equivalently here, S^a) and the Hubble parame- 
ter h. One can immediately notice in general a stronger 
degeneracy in all contour plots and, in particular, a more 
similar behavior of the bispectrum and power spectrum 
contours, especially for those involving Ha. This is the 
degeneracy that gets broken by acoustic features, as it is 
well known in the power spectrum case |66j. 

Focusing for instance on the Oa vs. h case in Fig. ^| 
and Fig. we can see that the bispectrum, by virtue of 
its several different triangular configurations, is remark- 
ably sensitive to features in the dark matter linear power 
spectrum such as the baryonic acoustic oscillations. The 
marginalized errors on individual parameters are, overall, 
larger when the BBKS transfer function is used. We no- 
tice, however, that the improvement provided by adding 
the bispectrum improves for parameters such as Ag, r, 61 
and the spectral index rig (about a factor of two better 
than the power spectrum alone) while it reduces for f2A- 



E. Dark Energy: wCDM models 



D. Baryon Acoustic Oscillations 

The same baryon acoustic oscillation features induced 
in the dark matter power spectrum |67l |6g| and re- 
cently seen in galaxy surveys l6J, |63| are expected to 
be present in the bispectrum [T^l and can also be used 
to help in determining cosmological parameters. Fig- 
ure IHl shows the ratio of the bispectrum to a featureless 
(no acoustic oscillation) bispectrum obtained from the 
BBKS fitting formula |4^ by setting the shape param- 
eter F = 0.175. Because the bispectrum scales as the 



We now extend the analysis performed above to in- 
clude the determination of the dark energy equation of 
state parameter w under the assumption of an homo- 
geneous dark energy component. We assume w to be 
constant. Introducing w, on which the growth function 
depends, leads to an increased degeneracy with the other 
parameters controlling the amplitude of the galaxy fluc- 
tuations such as As and 61 , which can be ameliorated by 
including bispectrum information. 

Table I Villi presents the expected errors on the vari- 



ous parameters for the two cases of k^ 



MS 



0.2/iMpc" 



18 



1.1 



fJTJTTTTTTTTTJTTTTTTTTTm 



1 -- 



0.9 
1.2 

1 

^ 0.8 

0.6 

1.2 
a" 
0.9 

1.2 

" 1 



i l|lllllllll|lllllllll ||j i| I m 



G 



0.8 


0.6 


0.8 


■^ 0.7 


0.6 


0.06 


= 0.05 



0.04 




] | lllllllll lllllllll 



il Ii 




ii 



lllllllll lllllllll 




II 



I Ii 





0.2 


XI 1 1 1 1 1 1 


1 1 1 1 II l_L 




0.1 


~ '7 


0~- — 







r 1 


/ ' - 




-0.1 


r '^ 


-_ 




-0.2 


-ri 1 1 1 1 1 1 


1 1 1 1 1 1 i-r 




0.9 1.2 
b, 



I lllllllll lllllllll [] i lllllllll lllllllll iii iiiiiiiiiiiiiiiiiiiiiiiiiiiii Di iiiii M ill iiiii : 




fStt 




III I I Il 



kn,ax = 0-3 h/Mpc 




IIIIIIIIIIIIIIIIIIII D: ||||| M ill lllli : 




1 1 1 1 



.1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1 1- 




: | I I I I I I I I I I I i : 




1 1 1 1 1 1 1 1 1 1 1 1 1 1 



P jrTTTTTTTTTTTTTTTTTTTTTTTTT 



Nir^ 



ii I Il 



0.6 0.7 0.80.9 1 1.10.6 0.8 1 1.2 0.9 1.2 0.60.8 1 1.2 0.6 0.7 0.8 
Ha "= K bi '^B h 



FIG. 10: Same as Fig. Elbut for a featureless (no baryon acoustic oscillations) transfer function. 



and 0.3/iMpc"\ We find that for fc^^f = 0.2/iMpc"^ 
the determination of the w parameter improves, over the 
W+P case, by 10% when the bispectrum is included, 
while by comparison a 20% if the LRG power spectrum 
is added instead. For fc^{f = 0.3 /iMpc^^ adding the bis- 
pectrum improves the determination of lu by 15%. These 
are mild improvements, but note that there are basically 
three different ways of getting below 10% errors on w 
by using the power spectrum of main and LRG samples, 
and the non-Gaussian information in the bispectrum and 



thus important consistency checks. Of course, the addi- 
tion of extra information such as type lA supernovae or 
weak gravitational lensing (apart from the latest CMB 
data) will tighten the constraints. 

In Fig. 111! we plot the marginalized likelihood func- 
tions for the case where the maximum wavenumber is 
k^S ^ o.3/iMpc"\ Note that unlike the ACDM mod- 
els previously considered, some of the maximum likeli- 
hood values of the first year WMAP data differ substan- 
tially from the chosen fiducial values for the LSS likeli- 



19 









TABLE VIII: Same as Table EH but for 


wCDM models 








W+P 


W+B 


W+P+B 


W+P+B (no mix. cov 


) W+P+B+Pi 


W+P+Pt 


W+Pl 




kiit = 


0.2/iMpc-i 








kO^f =0.2/iMpc-i 


Acjd 


0.0070 


0.0073 (-4%) 


0.0061 (15%) 


0.0058 (21%) ] 


0.0047 (49%) 


0.0050 (40%) 


0.0056 (25%) 


Aujt 


0.00109 


0.00113 (-4%) 


0.00094 (16%) 


0.00097 (12%) ] 


0.00084 (30%) 


0.00094 (16%) 


0.00102 (7%) 


AfiA 


0.0141 


0.0121 (17%) 


0.0102 (38%) 


0.0098 (44%) ] 


0.0076 (86%) 


0.0089 (58%) 


0.0101 (40%) 


Aris 


0.032 


0.032 (0%) 


0.026 (23%) 


0.027 (19%) ] 


0.022 (45%) 


0.025 (28%) 


0.029 (10%) 


AAs 


0.131 


0.129 (2%) 


0.107 (22%) 


0.112 (17%) ] 


0.092 (42%) 


0.109 (20%) 


0.119 (10%) 


Aw 


0.107 


0.103 (4%) 


0.098 (9%) 


0.096 (11%) ] 


0.084 (27%) 


0.089 (20%) 


0.095 (13%) 


At 


0.077 


0.075 (3%) 


0.067 (15%) 


0.068 (13%) ] 


0.058 (33%) 


0.066 (17%) 


0.071 (8%) 


Abi 


0.099 


0.132 (-25%) 


0.064 (55%) 


0.079 (25%) ] 


0.055 (80%) 


0.073 (36%) 


- 


A&2 


- 


0.072 


0.055 


0.063 ] 


0.051 


- 


- 


^f^LRG 


- 


- 


- 


0.101 


0.129 


0.14 


Ao-8 


0.076 


0.085 (-11%) 


0.049 (55%) 


0.057 (33%) ] 


0.042 (81%) 


0.056 (36%) 


0.066 (15%) 


Ah 


0.0183 


0.0164 (12%) 


0.0115 (59%) 


0.0123 (49%) ] 


0.0096 (91%) 


0.0119 (54%) 


0.014 (31%) 


AQt 


0.0028 


0.0025 (12%) 


0.0022 (27%) 


0.0022 (27%) ] 


0.0017 (65%) 


0.0019 (47%) 


0.0022 (27%) 




"'max — 


0.3/iMpc-^ 




fc,t2P =0.2;iMpc-i 


Aujd 


0.0068 


0.0056 (21%) 


0.0054 (26%) 


0.0053 (28%) ] 


0.0044 (54%) 


0.0050 (36%) 


Aujb 


0.00107 


0.00094 (14%) 


0.00087 (23%) 


0.00083 (29%) ] 


0.00080 (34%) 


0.00092 (16%) 


AQa 


0.0120 


0.0083 (45%) 


0.0082 (46%) 


0.0074 (62%) ] 


0.0066 (82%) 


0.0084 (43%) 


Aris 


0.031 


0.026 (19%) 


0.024 (29%) 


0.022 (41%) ] 


0.021 (48%) 


0.025 (24%) 
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0.096 (34%) 
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0.084 (54%) 


0.108 (19%) 
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0.091 (15%) 


0.085 (23%) ] 


0.080 (31%) 


0.088 (19%) 


At 


0.076 


0.066 (15%) 


0.061 (25%) 


0.059 (29%) ] 


0.053 (43%) 


0.065 (17%) 


Abi 


0.091 


0.100 (-9%) 


0.052 (75%) 


0.062 (47%) ] 


0.047 (94%) 


0.071 (28%) 


Afea 
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0.046 


0.042 


0.045 ] 


0.040 
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^f^LRG 


- 


- 
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0.085 


0.126 


Acrg 


0.070 


0.062 (13%) 


0.038 (84%) 


0.042 (67%) ] 


0.035 (100%) 


0.055 (27%) 


Ah 


0.0160 


0.0108 (48%) 


0.0091 (76%) 


0.0089 (80%) ] 


0.0083 (93%) 


0.0114 (40%) 


AQb 


0.0025 


0.0018 (39%) 


0.0018 (39%) 


0.0017 (47%) ] 


0.0016 (56%) 


0.0019 (32%) 



hood function. For instance, WMAP gives a marginal- 
ized probability distribution for w with a maximum close 
to —0.72, rather far from our fiducial value w — —1. 
This implies that part of the constraing power of the 
LSS statistics is spent in shifting the maximum to the 
w — —1. 

Finally, in Fig. El we show the contour plots for some 
selection of parameters. Comparing these results to the 
previous on ACDM models we see that the improvement 
brought by the bispectrum is increased in the case of 6i 
and as, and in fact the final error bars on these param- 
eters (and 62) are almost insensitive to including a more 
generic dark energy. This is good news as erg is one of the 
least known parameters and subject to tension between 
different data sets 



leters ana suDieci 

ilzllillzilzi 



On the other hand, the constraints on f2A, ^s, t and 
particularly Ud and Us are significantly worse than in the 
cosmological constant case. Regarding the behavior of 
the mixed covariance matrix, we see the same impact of 



beat coupling that we discussed before, the parameters 
responsible for the amplitude of galaxy fiuctuations As , 
T, as and bias parameters improve by the inclusion of the 
mixed covariance, though the behavior of the rest of the 
parameters is somewhat more complicated. Performing 
the same test as in Fig. |S1 in this case we see that ex- 
cluding the first fc-bin completely erases the effects due 
to beat coupling. 



VI. CONCLUSIONS 

We have provided a first detailed study on the infor- 
mation about cosmological parameters contained in the 
bispectrum of the galaxy distribution at large scales, pay- 
ing particular attention to the joint analysis with the 
power spectrum and their combination with CMB data. 
We have shown that the bispectrum has significant infor- 
mation on cosmological parameters and when combined 
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FIG. 11: mCDM models: marginalized likelihood functions for the six+four cosmological parameters assuming k„ 
0.3 /iMpc"^. For the bi parameter, the dotted line denotes the WMAP plus LRG power spectrum likelihood. 



with the power spectrum it is more complementary than 
combining power spectra of different samples of galax- 
ies, since non-Gaussianity provides a somewhat different 
direction in parameter space. Moreover, replacing the 
power spectrum with the bispectrum gives similar con- 
straints on cosmological parameters and can therefore 
serve as a consistency check. 

In order to properly combine bispectrum with power 
spectrum information, we worked out their covariance 



properties. Due to the effects of beat coupling [63|, the 
mixed terms in the covariance matrix are enhanced. We 
demonstrated that including this effect into the likelihood 
analysis provides a slight improvement on the error bars 
of cosmological parameters related to the amplitude of 
galaxy fluctuations. 

In the framework of flat cosmological models we 
showed that most of the improvement of adding bispec- 
trum information corresponds to parameters related to 
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FIG. 12; wCDM models: marginalized 95% contour plots for a selection of the cosmological parameters assuming fc„,ax 
0.3/iMpc-\ 



the amplitude and effective spectral index of perturba- 
tions, in particular f2„i (or Q\ = 1 — $!,„) and erg, which 
can be improved by factors of 1.5 to 2 and, interestingly, 
are presently among the least well known. In particu- 
lar, we showed that the uncertainties on trg are robust 
to relaxing the equation of state parameter w beyond a 
cosmological constant. This is good news as an is su bjec t 
to tension between different data sets 0, IZE IZ^ IZE UM ■ 
We also showed that the improvements are not directly a 



consequence of just constraining galaxy bias but of gen- 
uine information on cosmological parameters. 

As far as future theoretical work is concerned, the sin- 
gle most pressing issue is that of systematic errors in the 
predictions, which were addressed in (53 for a previous 
generation of galaxy surveys. Those methods, based on 
second-order Lagrangian perturbation theory (which we 
have used here) , are likely not enough given the expected 
statistical errors we derived in this work. Fortunately, 
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powerful methods based on first principles have recently 
become available |8Cl(. and together with numerical sim- 
ulations they should provide a sound theoretical basis. 
Work on this is in progress and will be reported soon. 
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APPENDIX: WMAP 3- YEAR UPDATE 

Shortly after the completion of the present work, the 3- 
year WMAP satellite observations, lij, became publicly 
available. For the reasons described below, the recent 
data is not simply an incremental improvement on the 
results presented in the main sections of the paper, and 
it is instructive to see the differences in the constraining 
power of the bispectrum. For this reason, and in order 
not to introduce substantial changes to preprint version 
of the paper, we add this appendix updating the results 
of sectionElto include the CMB likelihood corresponding 
to the new data. 

The improved analysis of the CMB polarization on the 
3-year data is responsible for significant differences with 
respect to the 1-year results in the errors on cosmological 
parameters as well as in their best values. Particularly 
significant to the present work is the much improved de- 
termination of the optical depth parameter r leading to 
a better constraints on the amplitude of primordial fluc- 
tuations As and on the spectral index n^.. 



TABLE IX: Fiducial values for the cosmological and bias 
parameters assumed for the LSS-I-WMAP 3-year likelihood 
analysis. 



Parameter 


ACDM models 


wCDM models 


LOd 


0.104 


0.105 


LJi 


0.0223 


0.0221 


^A 


0.765 


0.715 


Us 


0.95 


0.944 


As 


0.687 


0.715 


w 


-1 


-0.84 


T 


0.09 


0.09 


bl 


1 


1 


b2 








f^LRG 


2.17 


2.17 



We can generically expect that the smaller errors from 
the CMB analysis alone would make, for certain param- 
eters, the effect of adding other data sets less noticeable. 
In particular we can expect a reduced relative impact 
of including the bispectrum in the large-scale structure 
data analysis with respect to the parameters responsible 
for the amplitude of galaxy fluctuations. 

However, as we will show, the constraining power of 
the power spectrum and bispectrum joint analysis turns 
to other relevant parameters, particularly in the case of 
the wCDM models, where the uncertainty on the dark en- 
ergy equation of state introduces substantial degeneracies 
while the large-scale structure statistics are particularly 
sensitive to the late-time expansion via the growth factor 
D{p;a). 

Another reason to expect the improvement due to the 
bispectrum likelihood to be somehow smaller is related 
to the lower value of the amplitude parameter Ag and, 
consequently, of us • From Eq. Q and Eq. (|18|l one can 
derive the signal to noise for a given triangular configu- 
ration, assumed here equilateral for simplicity, so that, 
at large scales, we have 



(fc) 



B{k,k,k) 
AB{k,k,k) 



v/^W 



(52) 



The best value for the parameter A^ decreased from the 
1-year to 3-year WMAP analysis by almost 15%. On 
the other hand a lower value for amplitude of primor- 
dial fluctuations results as well in a lower value for the 
non-linear scale thereby making more robust our approx- 
imations based on perturbation theory. 

The definitions of the power spectrum and bispectrum 
likelihood functions, including their covariance proper- 
ties, are the same as the ones described in the previous 
sections. The sole exception that we will consider for 
this appendix regards the fiducial values assumed for the 
large-scale structure observables. To make the compari- 
son between power spectrum analysis and joint analysis 
as clear as possible, we take them to coincide with the 
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TABLE X: ACDM models: expected marginalized errors (68% CL) for WMAP3 (temperature and polarization, W) combined 
with the SDSS main sample power spectrum (P) and bispectrum (B) and with the LRG power spectrum (Pl). The percentage 
in parenthesis indicates the improvement over the analysis including the main sample power spectrum alone (W+P), numbers 
in bold indicate errors down by at least 1.5. In brackets we quote the W+P+B errors obtained by ignoring the mixed power 
spectrum - bispectrum covariance. Compare with Table IVTl for WMAPl. 





W+P 


W+B 


W+P+B 


W+P+B (no mix. 


cov.) W+P+B+Pl 


W+P+Pl 


W+Pl 




fc*lf = 


D.2/iMpc-^ 








kf^^ =0.2/iMpc-^ 


Auja 


0.0024 


0.0025 (-4%) 


0.0020 (20%) 


0.0020 (20%) ] 


0.0017 (41%) 


0.0018 (33%) 


0.0020 (20%) 


Aujb 


0.00062 


0.00063 (-2%) 


0.00058 (7%) 


0.00059 (5%) ] 


0.00056 (11%) 


0.00058 (7%) 


0.00060 (3%) 


A^A 


0.0088 


0.0077 (14%) 


0.0061 (44%) 


0.0060 (47%) ] 


0.0045 (96%) 


0.0052 (69%) 


0.0062 (42%) 


Aris 


0.0132 


0.0135 (^2%) 


0.0120 (10%) 


0.0123 (7%) ] 


0.0108 (22%) 


0.0116 (14%) 


0.0123 (7%) 


AA, 


0.038 


0.039 (-3%) 


0.034 (12%) 


0.036 (6%) ] 


0.034 (12%) 


0.037 (3%) 


0.037 (3%) 


Ar 


0.025 


0.025 (0%) 


0.023 (8%) 


0.024 (4%) ] 


0.023 (8%) 


0.025 (0%) 


0.049 (6%) 


A6i 


0.047 


0.059 (-20%) 


0.041 (15%) 


0.046 (2%) ] 


0.038 (24%) 


0.043 (9%) 


- 


Afea 


- 


0.067 


0.000 


0.060 ] 


0.044 


- 


- 


A&f«« 


- 


- 


- 


0.074 


0.081 


0.137 


Ao-8 


0.030 


0.031 (-3%) 


0.026 (15%) 


0.028 (7%) ] 


0.024 (25%) 


0.027 (11%) 


0.029 (3%) 


A/i 


0.0112 


0.0098 (14%) 


0.0086 (30%) 


0.0087 (29%) ] 


0.0073 (53%) 


0.0082 (37%) 


0.0089 (26%) 


Af^t 


0.00102 


0.00095 (7%) 


0.00080 (27%) 


0.00079 (29%) ] 


0.00068 (50%) 


0.00074 (38%) 


0.00082 (24%) 




fc*lf = 


].3/iMpc-^ 




kf^^ =0.2;iMpc-i 


Ac^d 


0.0023 


0.0020 (15%) 


0.0019 (21%) 


0.0018 (28%) ] 


0.0017 (35%) 


0.0018 (22%) 


AoJt 


0.00061 


0.00059 (3%) 


0.00056 (9%) 


0.00055 (11%) ] 


0.00054 (13%) 


0.00057 (7%) 


A^A 


0.0080 


0.0049 (63%) 


0.0049 (63%) 


0.0043 (86%) ] 


0.0039 (105%) 


0.0051 (57%) 


Aris 


0.0128 


0.0121 (6%) 


0.0111 (15%) 


0.0105 (22%) ] 


0.0102 (25%) 


0.0113 (13%) 


AA, 


0.038 


0.037 (3%) 


0.031 (22%) 


0.034 (12%) ] 


0.031 (22%) 


0.036 (6%) 


Ar 


0.025 


0.025 (0%) 


0.022 (14%) 


0.024 (4%) ] 


0.022 (14%) 


0.025 (0%) 


Afei 


0.046 


0.056 (-18%) 


0.038 (21%) 


0.042 (9%) ] 


0.036 (28%) 


0.042 (9%) 


A62 


- 


0.045 


0.036 


0.044 ] 


0.035 


- 


A&f«« 


- 


- 


- 


0.068 


0.080 


Ao-8 


0.029 


0.029 (0%) 


0.023 (26%) 


0.024 (21%) ] 


0.022 (32%) 


0.027 (7%) 


A/i 


0.0103 


0.0076 (35%) 


0.0071 (45%) 


0.0080 (65%) ] 


0.0065 (58%) 


0.0079 (30%) 


Af^t 


0.00096 


0.00073 (31%) 


0.00073 (31%) 


0.00068 (41%) ] 


0.00067 (43%) 


0.00074 (30%) 



maximum values of the WMAP 3-year likelihood. They 
are therefore different for ACDM and wCDM models and 
are reported in Table IIXI We make use of polynomial 
fits to the CMB likelihood functions determined in the 
same fashion as explained in section HV Bl but . this time, 
making use of the Monte Carlo Markov chains publicly 
available instead of evaluating the likelihood function di- 
rectly. 



1. ACDM models 

We present in Table Ixl the expected marginalized er- 
rors (1-<t) from the power spectrum and bispectrum like- 
lihood analysis combined with the WMAP 3-year data. 
As in section numbers in bold correspond to improve- 
ments larger than 50%, with the improvement factor de- 
fined in Eq. ifl^ . 



General considerations such as the dependence on the 
smallest scale (largest wavenumber k^^^) included are 
substantially unchanged. One can notice how in the 
fcmax = 0.2 /iMpc~^ case the improvement due to the bis- 
pectrum ranges from 7% for w;, to 44% for J^a, while 
in the fc„,^^ = 0.3/iMpc~^ case we have 9% to 63%. 
These results are just slightly lower than those obtained 
with the 1-year WMAP data. Instead, the improve- 
ments on the parameters responsible for the overall am- 
plitude of galaxy fluctuations and the spectral index are 
significantly reduced with respect to the analysis with 
WMAP 1-year data, this being a consequence, as men- 
tioned above, of the much better determination of the 
parameters t and As from the CMB data alone, whose 
error bars shrinked by a factor of 2 or more. This results, 
however, in much smaller expected errors on the bias pa- 
rameters, with the linear bias 61 determined to better 
than 4%. We can still conclude, that the combination 
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of CMB observations with the main sample power spec- 
trum and bispectruni obtains better constraints than the 
combination with power spectra from the two different 
samples considered in this work. 

Such considerations are reflected in the marginalized 
95% CL contour plots shown in Figure El where we can 
still observe the different directions of several degenera- 
cies of the power spectrum likelihood function with re- 
spect to the bispectrum one. 

We have checked that the conclusions derived from Ta- 
ble ^^ and discussed in section FV Bl hold as well for the 
WMAP 3-year data. In both tests, we can still observe 
better results in the W-|-P-|-B analysis, particularly for 
certain cosmological parameters such as ^m , again show- 
ing that the information provided by the bispectrum is 
not limited to the galaxy bias determination. 



2. Dark Energy: wCDM models 



Finally, in Table IXIl we present the expected marginal- 
ized errors (l-u) on the cosmological parameters for mod- 
els which allow for a homogeneous dark energy compo- 
nent with an equation of state parametrized by the con- 



stant w. 

As already mentioned, for these models the improve- 
ment due to the new CMB data set to the constraints on 
the parameters responsible for the amplitude of galaxy 
fluctuations is now less dramatic because of the degener- 
acy with w. In this case, therefore, we can expect a more 
relevant contribution of the bispectrum to the combined 
analysis because of the dependence of large-scale struc- 
ture statistics on the growth function D{p, a). 

We see indeed that, for fc,„ax = 0.3 /iMpc^^, includ- 
ing the bispectrum, the error on w improves by 46%, a 
significantly larger improvement compared to the previ- 
ous analysis using WMAPl (15%), or compared to the 
results obtained by adding the LRG power spectrum to 
the main sample one and WMAP3 (24%). We also find 
that in the most general case of W-t-P+B+P^ the same 
error gets better by a 54%. Similar observations hold as 
well for other parameters such as uJd and Qa while we 
notice less impressive results for Ag and n^. The corre- 
sponding two-parameter marginalized contour plots for 
the wCDM models updated to WMAP3 is shown in Fig- 
ure El where the difference in the improvements can be 
seen in more detail (compare to Fig ll2|) . 
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FIG. 13: Marginalized 95% contour plots for pairs of cosmological parameters in the ACDM models assuming k^^^ 
0.3/iMpc~^ and WMAP3 data. Compare with Fig.[7|for WMAPl. 
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TABLE XI: Same as TablelUbut for wCDM models. Compare with Table IVnTI for WMAPl. 





W+P 


W+B 


W+P+B 


IV+P+B (no mix. 


cov.) W+P+B+Pi 


W+P+Pl 


W+Pl 




fc*lf = 


3.2hMpc-^ 








A;f;ff =0.2/iMpc-i 


Awd 


0.0044 


0.0047 (-6%) 


0.0035 (26%) 


0.0036 (22%) ] 


0.0028 (57%) 


0.0032 (37%) 


0.0036 (22%) 


Aujb 


0.00060 


0.00061 (-2%) 


0.00058 (3%) 


0.00059 (2%) ] 


0.00057 (5%) 


0.00058 (34%) 


0.00059 (2%) 


AOa 


0.0104 


0.0092 (13%) 


0.0073 (42%) 


0.0072 (44%) ] 


0.0054 (93%) 


0.0062 (68%) 


0.0072 (44%) 


Aus 


0.0130 


0.0132 (-1%) 


0.0126 (3%) 


0.0126 (3%) ] 


0.0121 (7%) 


0.0124 (5%) 


0.0128 (2%) 


AAs 


0.037 


0.038 (-3%) 


0.035 (6%) 


0.036 (3%) ] 


0.035 (6%) 


0.036 (3%) 


0.037 (0%) 


Aw 


0.060 


0.058 (3%) 


0.045 (33%) 


0.048 (25%) ] 


0.040 (50%) 


0.046 (30%) 


0.051 (18%) 


At 


0.025 


0.025 (0%) 


0.025 (0%) 


0.025 (0%) ] 


0.025 (0%) 


0.025 (0%) 


0.025 (0%) 


Abi 


0.082 


0.101 (-19%) 


0.058 (41%) 


0.069 (19%) ] 


0.050 (64%) 


0.062 (32%) 


- 


Ab2 


- 


0.068 


0.051 


0.061 ] 


0.047 


- 


- 


Afef«« 


- 


- 


- 


0.091 


0.110 


0.121 


Ao-8 


0.048 


0.050 (-4%) 


0.034 (41%) 


0.038 (26%) ] 


0.029 (66%) 


0.036 (33%) 


0.041 (17%) 


Ah 


0.0141 


0.0123 (15%) 


0.0089 (58%) 


0.0096 (47%) ] 


0.0074 (91%) 


0.0093 (52%) 


0.0107 (32%) 


Afit 


0.00194 


0.00168 (15%) 


0.00134 (45%) 


0.00140 (39%) ] 


0.00116 (67%) 


0.00134 (45%) 


0.00150 (29%) 




ktii = 


3.3hMpc-^ 




kt^f =0.2/iMpc-^ 


Aud 


0.0042 


0.0038 (11%) 


0.0030 (40%) 


0.0029 (45%) ] 


0.0025 (68%) 


0.0032 (31%) 


Acjt 


0.00059 


0.00058 (2%) 


0.00058 (2%) 


0.00056 (5%) ] 


0.00056 (5%) 


0.00057 (4%) 


A^A 


0.0090 


0.0061 (48%) 


0.0059 (53%) 


0.0052 (73%) ] 


0.0048 (88%) 


0.0059 (53%) 


Aus 


0.0128 


0.0125 (2%) 


0.0124 (3%) 


0.0120 (7%) ] 


0.0120 (7%) 


0.0123 (41%) 


AAs 


0.037 


0.037 (0%) 


0.034 (9%) 


0.035 (6%) ] 


0.033 (12%) 


0.036 (3%) 


Aw 


0.057 


0.048 (19%) 


0.039 (46%) 


0.041 (39%) ] 


0.037 (54%) 


0.046 (24%) 


At 


0.025 


0.025 (0%) 


0.025 (0%) 


0.025 (0%) ] 


0.025 (0%) 


0.025 (0%) 


Abi 


0.077 


0.085 (-9%) 


0.049 (57%) 


0.057 (35%) ] 


0.043 (79%) 


0.060 (28%) 


Ab2 


- 


0.045 


0.039 


0.044 ] 


0.037 


- 


A&f«« 


- 


- 


- 


0.079 


0.108 


Ao-8 


0.046 


0.040 (15%) 


0.027 (70%) 


0.029 (59%) ] 


0.024 (92%) 


0.035 (31%) 


Ah 


0.0124 


0.0084 (48%) 


0.0070 (77%) 


0.0069 (80%) ] 


0.0063 (97%) 


0.0089 (39%) 


AQb 


0.00173 


0.00122 (42%) 


0.00114 (52%) 


0.00113 (53%) ] 


0.00106 (63%) 


0.00130 (33%) 
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FIG. 14: uiCDM models: marginalized 95% contour plots for a selection of the cosmological parameters assuming fc„,ax 
0.3/iMpc-\ Compare with Fig.IHlfor WMAPl. 
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